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Abstract 

Discrete Cosine Transform (DCT) is very important in image compression. Classical 1-D DCT 
and 2-D DCT has time complexity 0(N log 2 N) and 0(N 2 log 2 N) respectively. This paper presents 
a quantum DCT iteration, and constructs a quantum 1-D and 2-D DCT algorithm for image 
compression by using the iteration. The presented 1-D and 2-D DCT has time complexity 0(y/N) 
and O(N) respectively. In addition, the method presented in this paper generalizes the famous 
Grover's algorithm to solve complex unstructured search problem. 
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I. INTRODUCTION 



A. The Conception of Digital Image 



The result of ^sampling and quantization of a monochromic N-by-N image is a matrix of 
real numbers [3, 
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/(JV-1)0 f(N-l)X ••• f(N-X)(N-X) . 



(1) 



Each element of this matrix is called pixel. The pixel fy is a gray-level value. The value 
ffi is proportional to brightness or energy. 



Let column vector 



fi — {foil flh ■••■> f(N-l)-i) J 

, where T represents the transpose of a vector. 
The digital image F can be expressed as formula 



fo fi ■■■ /at-i 
Fig. 1 illustrates the conception of digital image. 

Insert Figl Here 



(2) 



(3) 



B. Introduction of Classical Discrete Cosine Transform (DCT) 

1. One-Dimensional Discrete Cosine Transform (1-D DCT) 

The one-dimensional discrete cosine transformQ (1-D DCT) and inverse DCT are 
defined as 

N-x 

Q-uJn COS 2J\r 

ra=0 

j? _ V> (2n+l)u7r (4) 

In / j ^u^u COS 



where 



n=0 

(0 < u,n <N- I) 
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^/N if 1 < u < N 
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The NxN DCT matrix is 
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<xn-i cos ^ 2jV ; • ■ • a N _i cos i y 



Let vector 



(5) 



~> 1.147T 3.M7T 

D u = a u (cos -^r, COS -^r, COS 

, where < u < N 

The NxN DCT matrix can be expressed as 



(2iV-l).wr, 
2N ' 



(6) 



(2n + l)wr 
D = (a u cos — ) NxN 
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Equation (4) can be expressed as 
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where "•" expresses inner product between two vectors. 



(8) 



2. Two-Dimensional Discrete Cosine Transform (2-D DCT) 

The definition of the two-dimensional DCT 0, H (2-D DCT ) for input image F and 
coefficients matrix C (output) is 



N-XN-l 

E\—^ n 7r(2m+l)p 7r(2n+l)g 

/ j Ivan COS 2 JV COS 2 JV 
m=0 n=0 



r _ 7r(2m+l)p 7r(2n+l)(j l y J 

Jmn 2 / 2 / O^p^qCpq COS 27V COS 27V 
m=0 n=0 

(0 < p,q,m,n < N) 

2-D DCT is closely related to the discrete Fourier transform. It is a separable linear 
transformation. That is, the result of 2-D DCT may be obtained by first taking transforms 
along the columns of F and then along the rows of that result 0, 0|- That is [2,0], 

C = DFD (10) 
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In image compression, the input image is divided into 8-by-8 or 16-by-16 blocks, and the 
2-D DCT is computed for each block. The DCT coefficients are then quantized, coded, and 
transmitted. The receiver decodes the quantized DCT coefficients, computes the inverse 2-D 
DCT of each block, and then puts the blocks back together into a single image. 

Almost all of digital films such as VCD and digital pictures such as JPEG files are 
compressed by using DCT currently. Real-time compressing and high compression ratio are 
the main research topics of image compression [2] . The time complexity of classical 1-D DCT 
and 2-D DCT is 0(JVlog 2 iV) and 0(iV 2 log 2 iV) respectively 00. The more N is big, the 
more the compression ratio is high if the whole input data sequence is Markov chain and the 
order of chain is bigger than N, but the running time increases drastically so that real-time 
compressing is impossible. That is, finding a fast algorithm for large N is significant. 

3. Two Properties of DCT of Image 

• Property 1: 2-DCT of image is an energy conservation transform 0,0, i.e., 

N-1N-1 N-lN-1 
m=0 n=0 p=0 q=0 

1-D DCT is same too. 

The property is utilized as halt criterion of algorithm in this paper. If we find some DCT 

\ i \ i 

coefficients (e.g., c 00 , Coi, c w ) such that J2 J2 (f mn ) 2 ~ (c o) 2 + (c i) 2 + (ci ) 2 , we will halt 

m=0 n=0 

algorithm. That is, coo,coi,cio keep the information of image approximately. 

• Property 2: For typical images, many of the DCT coefficients have values close to 
zero; these coefficients can be discarded without seriously affecting the quality of the 
reconstructed image 0, . 

Fig. 2, Fig. 3 and Table 1 illustrate the property 2 [3]. The property is utilized to design 
quantum algorithm in this paper. 

Fig2 Fig3 

(suggested minification: 60%) (suggested minification: 60%) 

Table 1 



C. The Classical Parallel Circuit for Computing Inner Product 
Notation 1 tA denotes the running time of addition a + b or subtration a — b. 
Notation 2 tu denotes the running time of multiplication a x b. 

Notation 3 to denotes the running time of division ~. In general, the time of division is 
little bigger than the time of multiplication. 

Notation 4 tc denotes the running time of comparison a < b or a > b. 
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In classical computer, there exist constants mi and such that tu > m\tA and to > 
rri2tM- The values of ni\ and ni2 relate with the material architecture of central processing 
unit (CPU). Different CPU has different values of mi and 777.2 0|. For the convenience of 
computation, without loss of generality, we suppose 

t M > 25t A , t M = 2t M and t A > £ c 

The assumption does not affect the time complexity because time compelxity is only 
associated with the order of parameter N. E.g., 0{\y/N) = 0{VN) and 0(| \/N) ^ 0{N). 
Let N — demensional vector 

x = (x , X\, ■ ■ ■ , xn-i) t 

y = (yo, yi, ■■■ ,yN-i) T 

The inner product is defined as 

JV-l 

x my = ^Xk x y k 

k=0 

The classical parallel circuit for computing inner product can be designed easily. Fig. 4 
illustrates a realization of the classical parallel circuit. 

Insert Fig4 Here 

Notation 5 tj denotes the running time of parallel computation of inner product. 

Clearly, ti = 1 X t M + Rog 2 N]t A 

when N < 2 4x25 , t / <lxi M +4x (25t A ) < 5t M 

Because = 2 100 is a very giant number, we regard £/ as 

tj < 5t M (11) 
D. Introduction of Quantum Computation 

Recent research shows that realizing image compression by using quantum computer is 
possilbe. For example, Pang recently presents a quantum algorithm to realize the encoding 
of Vector Quantization, which is very faster than classical encoding algorithm 0. This 
paper presents a qunatum algorithm to realize DCT for image compression. 

1. Introduction of Quantum Bit (qubit) 

Just as classical bit has state - either or 1 - a qubit also has a state. Two possible 
states for a qubit are the states |0) and |1), which as you might guess correspond to the 
states and 1 for a classical bit. Notation like '| )' is called the Dirac notation, and we'll be 
seeing it often, as it's the standard notation for states in quantum mechanics. The difference 
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between bits and qubit can be in a state other than |0) or |1). It is also possible to form 
linear combinations of states, often called superpositions 

|V) = a|0>+/3|l> 

The numbers a and j3 are complex numbers. Put another way, the state of a qubit is a 
vector in a two-dimensional complex vector space. The special states |0) and |1) are known 
as computational basis states, and form an orthonormal basis for this vector space. 

We can examine a bit to determine whether it is in the state or 1 in classical computer. 
By contrast, when we measure a qubit, we get either the result 0, with probability \a\ 2 , or 
the result 1, with probability 2 . Naturally, |a| 2 + |/3| 2 = l. In general a qubit is a unit 
vector in a two-dimensional complex vector space jfj. 

We explain the superpositions = a\0) + /3|1) by analogy with sonic wave as following. 

Suppose there are three persons Alice, Bob and you in a closed room. Alice and Bob 
speak in a sample wave \A) = f A {t) = e muJt and \B) = /#(£) = e muJt respectively, where 
m ^ n and they are both integers. We can distinguish Alice from Bob because the two 

r — 

sample waves are orthogonal (i.e., (A\B) = J " e l{ - n - m > ut dt = 0). When Alice speak in the 
closed room, your ears will receive a sonic wave g A (t) = I A e^ A \A), where I a is the amplitude 
of the wave and the phase <p A is cause by the distance between Alice and you. If Alice and 
Bob speak simultaneously, your ears will receive a superposition: 

\iP AB ) = I A e* A \A) + Ise^B) 

Let 

I = VUa) 2 + (Ib) 2 
a A = 1 -fe^ A and (3 B = 

Thus, 

\i/>ab) = I(<*a\A) + p B \B)) 

Your ears can distinguish Alice's voice from the superposition \iPab)- That is, '+' implies 
two sonic wave \A) and \B) exist in the superposition simultaneously and they can be 
distinguished from the superposition. If Alice speaks very aloud (i.e., \a A \ 2 — > 1 or \/3b\ 2 — > 
0) , you will always hear Alice's voice. This case is analogous with the case \a\ 2 1 of 
quantum computation. If \a\ 2 — > 1, you will get the result always, with probability | ct | 2 ~ 
1. This property is utilized to design quantum algorithm such as Grover's algorithm. You 
can operate the two distinguished sample wave \A) and \B) simultaneously. For example, 
you can send the voice \4>ab) i n a radio and change Alice's volumes and Bob's volumes 
simultaneously by pushing the volume button on radio. That is, performing once operation 
causes the changing of two sonic waves simultaneously. This case is analogous with quantum 
parallelism. Fig. 5 shows the analogies between quantum superpositions and sonic wave. 

Insert Fig5 Here 



2. Operation of Computation Acting on Qubits 

Classical computer circuits consist of wires and logic gates. The wires are used to carry 
information around the circuit, while the logic gates perform manipulations of the informa- 
tion, converting it from one form to another. It is the fundamental of classical computation 



5 



that classical computer circuits can realize the operations of Boolean algebra. For example, 
classical NOT gate makes and 1 states interchanged. The operations of Boolean algebra 
can also be realized on quantum computer by utilizing single qubit gates and controlled- 
NOT gates. For example, quantum NOT gate takes the state = a\0) + /3\1) to the 
corresponding state in which the role of |0) and |1) have interchanged. All digital operation 
can be realized by utilizing unitary operation that is the combination of some single qubit 
gates and controlled-NOT gates p. 



3. Quantum Parallelism 

Quantum parallelism allows quantum computers to evaluate a function f(x) for many 
different values of x simultaneously. The power of quantum computation is due to the fact 
that the state of a quantum computer can be a superposition of basis states and we can 
perform an operation on multiple quantum states simultaneously. For example, suppose 
f(x) : {0,1} — > {0,1} is a function with a one-bit domain and range. We need at least 
two times calculating for obtain the values f(0) and f(l) on classical computer. For arbi- 
trary function f(x), there is quantum circuit Uf that can transform \0,y) and \l,y) into 
|0, y® /(0)) and \l,y © /(l)) respectively by performing only one time calculating, where 

© indicates addition modulo 2. That is, |0)l ^ 1)l;/) , w here '+' im- 

plies two states |0) \y) and |1) \y) exist in the superposition of states simultaneously. The 
formula implies also that the two states |0) \y) and |1) \y) are converted to |0) \y®f(0)) and 
11) Jyffi/(1)) simultaneously. That is, the values /(0) and f(l) are calculated simultaneously 
00. 



4- Introduction of quantum Fourier transform (QFT), Klappenecker's DCT and Latorre's 
Quantum Representation of Image 

The quantum Fourier transform (QFT 0,[t|) on an orthonormal basis |0), |1), • • • , \N— 1) 
is defined to be a linear operator with the following action on the basis states, 



1 N ^ 

QFT\j) = J2 e2mjk/N \ k ) 



N-l 

Equivalently, the action on an arbitrary state may be written 

N-l N-l 

QFT(J2^\j)) = $>,(QFT|j» 

j=0 j=0 

Similar to other unitary operation, QFT is unitary operation that only acts on basis 
states. It's an error opinion that unitary operation can act on coefficients of basis states. 
Figure 6 illustrates that QFT only acts on basis states |0) and |1). 

Insert Fig6 Here 

Klappenecker presents a method to realize DCT of types I, II, III, and IV on a quantum 
computer by utilizing QFT [8[. 
Define DCT of types I as 
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C n '■= \l (TF)N c os^r]i,i=o...iv 



2 

The discrete sine transforms (DST) of type^s I denoted by S* N is denned accordingly, 
Let discrete Fourier transform (DFT) be 



N-l 



N 



2nijk/N 



)]k,l=0.. 
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Klappenecker's DCT derives from QFT and depends on QFT. Indeed, the DST and DCT 
can be recovered from the DFT by a base change 

T£ . F 2 n ■ T N = C T N © iS N 

Since efficient quantum circuit for the DFT (i.e., QFT) are known, it remains to find an 
efficient implementation of the matrix T/y. A quantum circuit is proposed by Klappenecker 
to realize the matrix T/y. This is the primitive idea of Klappenecker's DCT 

The result of QFT or Klappenecker's DCT seems to indicate that quantum computer can 
be used to very quickly compute the Fourier transform, which would be fantastically useful 
in a wide range of applications. However, that is not exactly the case; the Fourier transform 
or Klappenecker's DCT is being performed on the information 'hidden' in the amplitudes of 
the quantum state. This information is not directly accessible to measurement. The catch, 
of course, is that if the output state is measured, it will collapse each qubit into the state |0) 
or |1), preventing us from learning the transform result directlyjy]. In addition, the contents 
of section 2.1 in this paper shows that Klappenecker's DCT cannot be applied to DCT 
of image compression too. Klappenecker's DCT is useful on many other signal processing 
maybe. 

Latorre presentes a novel quantum representation for image compression The en- 
tropy of images is, in general, very large. An image with large entropy is hard to compress. 
The idea in Latorre's paper is to keep the largest eigenvalues of the Schmidt decompositions 
when the picture is written in a renormalization group manner, that is, the largest contribu- 
tion to the entropy in that basis. Latorre's algorithm is nice but it is not competitive with 
jpeg. The reason is that jpeg uses details of how human see. The quantization table used 
in jpeg is tailored to the human eye. Unless there are quantum quantization methods are 
incorporated in Latorre's algorithm, there is no way it is as efficient as jpeg. 
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E. Introduction of Grover's Algorithm and Boyer's Algorithm 



The algorithm by Grover solves the problem of searching for an element with a unique 
index xq in a list of N unsorted elements, similar to searching a database like a telephone 
directory when we know the number but not the person's name [jj. Grover's quantum 
searching algorithm [j| takes advantage of quantum mechanical properties to perform the 
search with an efficiency of order 0{\fN) p. 

To implement the quantum search we need to construct a unitary operation that discrim- 
inates between the marked item xq and the rest. The following function, 



/(*) 



if x 7^ xo 

1 if X = Xq 



and its corresponding unitary operation (i.e., black box or oracle), 

0\x) = 

will do the job. 

The Grover iteration is defined as 



G = (2\V)(V\-I)0 (12) 



N-l 

on N — dimensional Hilbert space, where \^f) = -4= \x). 



x=0 



The Grover iteration can be regarded as a rotation in the two-dimensional space spanned 
by l^) and \x ). To see this it is useful to let \a) = ^—^ Yl \ x ) indicate a sum over all x 

which are not solutions to the search problem. Thus, 



l*> = V^ |a) + V^ |xo) 

so the initial state of the quantum computer is in the space spanned by \xq) and \a). Let 



sin 0/2 = y/T/N 

As Figure 7 shows, G is a rotation in the two-dimensional space span{\xo), |ct)}, rotating 
the state by 6 radians per application of G. Repeated application of the Grover iteration 
rotates the state vector close to |a;o). When this is occurs, an observation in the computa- 
tional basis produces with high probability p. This is the key of Grover's algorithm. 

Insert Fig7 Here 

Eli Biham's papers show that Grover's algorithm is suitable for the case of an arbitrary 



initial amplitude distribution |lOt lllj. Gui-Lu Long improves Grover's algorithm to 100% 
successful probability even the number is very small [l^, 3 • 

Michel Boyer, et., al. present the following algorithm for the case that the number of 

solutions t is unknown 

Algorithm 6 Boyer's algorithm 



S 



1. Initialize m = 1 and set A = 6/5. 

(Any value of A strictly between 1 and 4/3 would do.) 

2. Choose j uniformly at random among the nonnegative integers smaller than m. 

3. Apply j iterations of Grover's algorithm starting from initial state |^o) — Yli 

4. Observe the register: let i be the outcome. 

5. If T[i] = x, the problem is solved: exit. 

6. Otherwise, set m to min(\m, x/iV) and go back to step 2. 

Boyer's algorithm is derived from Grover's algorithm and has time complexity 0(\J^-). 

The case of t = is handled by it. It is a very smart algorithm becuase it saves many 
quantum circuits comparing wiht the other algorithms by using quantum counting technique 
that will needed exponential order circuits p. Boyer's algorithm is utilized in this paper. 
The one of task of this paper is to design appropriate quantum iteration (named quantum 
DCT iteration in this paper) according to DCT properties to substitute for the Grover 
iteration in Boyer's algorithm. 

In this paper, Each ket represents a register rather than a single qubit. In this paper, 
similar to other quantum algorithm 0, 0] , ancilla qubits are ignored. 



II. THE REPRESENTATION OF IMAGE BY USING QUANTUM STATES 
A. Classical Data Structure of 1-D DCT 

For a given vector x = (xq, xl, ■ ■ • , x^v-i), we can declare a BYTE array "BYTE x[N}" 
to store it by using c language, where c language is compiler language of classical computer 
0: 

x[0] = Xo, x[l] = x 17 x[N — 1] = xn-i 
There is a logical mapping to associate subscript with component of vector x: 

Mapping : i i — ► x[i] ^ „s 

(0 < i < N) ( 16) 

The logical mapping is necessary because it associate data with the corresponding logical 
address. CPU accesses value x[i] according to the subscript i (i.e., logical address). 

The mapping is done by memory- management unit (MMU) of Operating System (l6| . 
The operation of access data is a very very fast operation so that the time of access can be 
ignored when designing algorithm. 

Fig8 illustrates the logical mapping. Fig9 illustrates the physical realization of the logical 
mapping 

Insert Fig8 Here 



Insert Fig9 Here 

Similar to vector x, the vectors (equation 2), c (equation 8), D u (equation 6), matrix D 
(equation 7) and matrix F (equation 1) can be stored in array respectively, and the Operating 
System of classical computer will establish the mapping (equation 13) automatically }l6^ . 



9 



For example, we declared a two dimensional array "BYTE array Image[N][N}" to store 
matrix F. The mapping between position (i,j) and pixel value is defined as 

Mapping : i— > array/ma<7e[i][7'] = (14) 

The above mapping (equations 13 and 14) should be also kept in quantum computation so 
that arbitrary component of vector or matrix is associated with the corresponding subscript. 

By the definition of DCT (equations (4) and (9)), QFT and Klappenecker's DCT cannot 
both keep the mapping. Therefore, More suitable quantum data structure is required in order 
to keep the mapping and harness the power of quantum computation for image compression. 



B. The Quantum Representation of Image 

1. Data Structure of Quantum Representation of Image 

To keep the mapping in equations (13) and (14), the following database technique is 
presented to represent image F (equation (1)) in this paper: 

First, all vectors /j = (foi, fu, f(N-i)i) T (0 < i < N) are stored in a database. Each 
vector is regard as a record with unique index i. 

Second, all vectors are loaded into CPU simultaneously and form the superposition of 

JV-l -> 

states -j=\ancillal) ( ^2 \i)\fi))\ancilla2) by using quantum addressing scheme and unitary 

i=0 

operation LOAD. 

LOAD operation that is denoted by Ul is defined as 

|z>|o> - - - |o> ^ |z)|oe/«)---|o 

, where ffi denotes addition modulo 2, that can be realized by utilizing controlled NOT 
operation 

In vector notation, 

|i>|o)&|i)|oe£) (is) 



LOAD operation is the basic operation of quantum computer (p, chapter 6). 
Figure 10 illustrates the representation of image by using quantum states. 

Insert Fig 10 Here 

The proposed representation of image in this paper keep the mapping in equa- 
tion (14) so that subscript (J, i) is associated with corresponding pixel value fa. State 

N-l -» 

■^=\ancillal) ( \i)\fi))\ancilla2) is entanglement state when ancillal and ancillal are con- 

i=0 

stants. Therefore, if we obtain value i from the second register, we will get the unique 
mapping vector /j in third register. Thus, the mapping is kept. 
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2. The Time of Load Operation Ul 



Notation 7 ti denotes the time of performing one time operation \i)\0) i— > |i}|0 © fji), 
where fji is a component of vector fi . 

Notation 8 Time(U) denotes the time of performing one time operation U. 

The time of LOAD operation Ul in equation (15) is denoted by Time(L r i). 

The time of data access (i.e., loading data into registers from memory) is ignored in 
classical algorithm when design algorithm. It is clear that the most efficient possible al- 
gorithm is in this model of computation chapter 6, section 6.5 or |l6T|). Operation 
|z)|0) I— > |z)|0 © fji) is CNOT operation, and it is a very very fast operation so that it's 
running time can be ignored when analyze the time complexity of algorithm ( 6], chapter 
6). Because each component fji is saved in memory independently, each component can be 
load into register simultaneously by utilizing parallel circuit in classical case. Furthermore, 
the classical circuit can be translated into a quantum reversible circuit. j6[ Therefore, in 
general, 

Time{U L ) =t L ~0 (16) 

III. QUANTUM 1-D DCT ITERATION 

A. The Design of Oracle O inner , 0' inner and O f 

All of vectors Di in equations (5) and (6) are stored in a database (FiglO). Each vector 

Di has a unique index i, where < i < N . 

An arbitrary vectors Di is loaded into registers according to the corresponding index % 
by using unitary operation Ul (FiglO): 

|«>|/3)|z)|/)|0)|0)|0)^|a>|/3)|z)|/)|A)|0)|0) (17) 
where a and (3 is two input parameters. 



1. Oracle O inner and 0' inner 

It's easy to construct a classical parallel circuit to calculate inner product (Section I-C 
and Fig. 4). Furthermore, the classical circuit can be translated into a quantum reversible 
circuit p, that has the same time complexity with it's corresponding classical circuit. 

We design oracles Oi nner and 0' inner to compute the inner product between two vectors Di 

and /. 

|a> \(3) |i> |?) | A) |0> |0> \a) \(3) |<> |/)|A)|0) I (A • If) (18) 

l«> 1/5) \i) If) I A> |0) |0) °*r \a) \(3) \i) |?) | A) | A • /> |0> (19) 
Fiure 11 illustrates the oracles Oi„ VPr - 
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Insert Fig 11 Here 

By equation (11), oracles O inner has time complexity 

) = t T + 1 x t M < 6t M (20) 

Oracles 0' inner has time complexity 

Time{0' mner ) =tj< 5t M (21) 

2. Oracle O f 

Define a function / as following 



lifa<(D t *f) 2 <[3 
otherwise 



We design other oracle denoted by Of that is used to mark the retained DCT coefficients 
(Table 1): 



|a>|/3)|i)|/)|A)|0)|(A.?) 2 )^ 

(-l)/W|a)|/3)|i)|/)|A)|0)|(A./) 2 ) (22) 

Figure 12 illustrates the oracles Of. 

Insert Fig 12 Here 

Oracle Of has time complexity . 

Time{O f ) = 2t c (23) 

B. Quantum 1-D DCT Iteration G DC t 

Let (Oinner) -1 denote the inverse transform of O inn er- 
Let (Ul)' 1 denote the inverse transform of U L . 

{OinnerU]^) (J-^l) (Oinner) 

Definition 9 By equations (17), (18) and (22), quantum 1-D DCT iteration G DC t is de- 
fined as 

G D CT = (2|0<£l " I){U L )-\O m ner)- l OfO m nerU L (24) 

N-l 

, where |f> = 7)7 £ 10- 

i=0 

12 



We have 

(U L )-\O inner )- l O f O mner U L \a) \(3) \i) | /) |0> |0) |0) = 

(-lp)|a)|/3)|i)|/)|0)|0)|0) (25) 

Equation (25) shows that the state |a)|/3)|«)|/)|0)|0)|0) is the eigenstate of unitary oper- 
ation (UL)~ l (Oi nner )~ l O fOi nner UL and the eigenvalue is (— This is a key of applying 
two or more oracles to perform complex search. 

Figure 13 shows the action of unitary (UL)~ 1 (O inner )~ 1 OfO inner UL. 

Insert Figl3 Here 

Let span{\i) <g> \vector) <g> \inner product)} be the global space that spanned by N — 
dimensional vector \vector), and the value of inner product \inner product) . 
Let span{\i) | < i < N} be the subspace that is the span of all states \i). 

Theorem 10 (Figure 14, 15) Gdct is a rotation on N — dimensional subspace span{\i) \ 

N-l -> 

< % < N}. For initial state |^o) = ~m E (l")!/ 5 )!^)!/)! )! )! )); the rotation angle is 



6 = 2 arcsin 




where t is the number of solutions. 



Proof 11 For arbitrary input state 

N-l 

\inpnt) = j2M\a)\m)\j)\m\o)) 

i=0 

where Ai is the probability amplitude of the state and 

N-l 

E i A - 1 2 = 1 

By equation (25), we have 
G DC T\input) = | a) 

N-l 

{(2|O(e|-/)(E(- 1 ) /( ^^)}l?)l )l°)l ) ( 26 ) 

i=0 



N-l 

Equation (26) shows that unitary operation 2|£)(£| — / only acts on state E ( — iK^^K) ■ 
Suppose that the solutions set is 

S = {\a)\m\l)\0)\0)\0)\a<(D t .Jr<(3} 
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All solutions form a superposition 



i r > = 4^ E i M ) = E i*» i/>io)|o>io) 



S \\u)es v\ S \ . s 



Lei 



|a) = -=J=y>) = 



1 



la if 

j-s 



T/ie initial state \ip ) is rotated in the subspace span{\a) ,|r)} fry towards the superposi- 
tion |t) o/ a// solutions to the search. 

Initially, it is inclined at angle | from \o~), a state orthogonal to |r). 

5?/ equation (25), the product operation ( y UL)~ 1 (O inner )~ 1 OfO inner UL reflects the state 
about the state \a). Then the operation 2|£)(£| — / reflect it about 

Notice that |r) G span{|i)} and \a) G span{|i)}. 

Therefore, similar to Grover's iteration, the quantum 1-D DCT iteration Gdct is a 
rotation on subspace span{\i)}, and the rotation angle is 9 = 2arcsin 

Figl4 illustrates the action of a single Gdct- 

Figl5 shows that Gdct is equivalent to a rotation on subspace span{\i)}. 

Insert Fig 14 Here 



Insert Fig 15 Here 

l-D DCT iteration Gdct has time complexity 

Time(G dct) = 2Time(U L ) + 2Time{O inner ) + Time(O f ) 
By equations (16), (20) and (23), when N < 2 100 , 

Time(G D cr) < 2t c + l2t M 
Because t M > 25t A 3> t c , we have 

Time(G dct) < 13t M (27) 

The theorem shows that the quantum 1-D DCT iteration Gdct and the Grover iteration 
G = (2|^) (ty \— 1)0 have the same property that they are both a rotation on space. However, 
comparing with Grover iteration that has only one oracle, Gdct comprises two oracles Oi nner 
and Of so that it can perform complex search. Equation (25) is a key of applying two or 
more oracles to perform complex search. If there is no equation (25), applying two or more 
oracles is no significant. 
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IV. QUANTUM 1-D DCT OF IMAGE COMPRESSION 

A. Subroutine 1: Finding a DCT Coefficient Cj G {cj | a < (cj) 2 < f3) 

Subroutine 1 is similar to Boyer's algorithm, and the Grover iteration in Boyer's algorithm 
is substituted by the iteration Gdct- 

Algorithm 12 Subroutine 1 

Input parameters: a, (3 

Outputs: subscript z'o and coefficient Cj ; a global Boolean variable IsSolution. 
Stepl. Initialize m — 1 and A = |; IsSolution = FALSE; 

Step2. Choose j uniformly at random among the nonnegative integers not bigger than 

m. 

Step3. Generate the initial state 



1 N ~ l 
* > = — |a>|/3>(^|z>)|/>|0>|0>|0> 



JV-l 

This can be achieved in 0(log2N) steps using a [log 2 N] — bit Hadamard transfor- 
mation. 

Step4. Apply j iterations of Gdct starting from the state \^o). i.e., 

|*i> = (G D ct) j \^o) 

Step5. Perform computation . 

Step6. Observe the third and the fifth register: let i , c io be the outcome respectively. 
Step7. If c io G | a < (q) 2 < j3}, Let IsSolution = TRUE, and exit; otherwise, go 
to step 8. 

Step8. set m to min{Am, \/N} and go to step 2. 
Figure 16 shows the processing of subroutine 1. 

Insert Fig 16 Here (suggested minification: 78%) 



B. The Correctness of Subroutine 1 



When 2 s < N < 2 1U0 , | - 1| < 6.5 x 10" 4 , 99.6% < 1 - f < 1. That is, when 

2 8 < iV < 2 100 , Grover's algorithm is still an efficient possible algorithm and the success 
probability p = 1 - o(^) > 99.6% 0. Long improves Grover's algorithm so that the 
improved algorithm with 100% success probability even iV is very small [Hill. 

The theorem in this paper demonstrates that Gdct is similar to the Grover iteration. 

Therefore, subroutine 1 is correct. 
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C. The Time Complexity of Subroutine 1 

Let k presents the iteration steps of applying Gdct- 
Subroutine 1 has time complexity 

Time(subroutinel) < K,[Time(G dct) + Time{0' inner ) 

+ |~log 2 N]Time(Hadamard) + 2t m + 3t c ] 

, where Time(Hadamard) represents the time of Hadamard transform acting on a single 
qubit. 

when N < 2 100 , we have 

log 2 NTime(Hadamard) xs 
Because tu > 25t^ ^> t c , by equations (21) and (27) 

Time(subroutinel) < 21k£m 
Subroutine 1 is similar to Boyer's algorithm. Hence 

Therefore, 

[N 

Time(subroutinel) = 0(y — )t M (28) 

D. Quantum 1-D DCT of Image Compression 

Algorithm 13 Quantum 1-D DCT of Image Compression (QDCT1) 

Input: {-DJO < % < N}: The line vectors of DCT matrix (equation 5, 6, and 7). These 
line vectors are stored in database with indices i (FiglO); 

Vector / (equation 2) 

Parameters:A_E, a, (3, threshold value e > 0; 

nCounti (0 < i < N): Counting the repeated times of coefficient q that is found when 
repeatedly applying subroutine 1. 

nM ax Repetition: The maximum number of allowed repeated times of applying subrou- 
tine 1. 

nS: The number of solutions that have been obtained. 

Output: The big DCT coefficients and corresponding indices (i ,c io ), (i 1: c ii: ), 

R 

(iR,c iD ) such that fc=0 < e, where 

coefficients are set to zero (Table 1). 

Step (Initialize parameters): AE 

0, where < i < N. 



/II 2 = (/•/)= E(/i x Si)- The other 

i=0 

|]7|| 2 , a = ^.,P = AE; nS = 0; nCounti = 
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The proposed parallel circuit in section I-C (figure 4) can be applied to calculate the 

value ||7|| 2 = (?•?)• 
while(^- > e) 

ll/ll 2 

{ 

Step 1: Apply subroutine 1 to find a coefficient (ik,c ik ); 
Step 2: If IsSolution = TRUE and nCounU = 0, 

then AE = AE - c ik x c ik ; nS = nS + 1; a = J^ s , (3 = AE; nCount ik = 
nCounti k + 1; 

If IsSolution = TRUE and < nCounti k < nM ax Repetition, continue; 

If IsSolution = TRUE, and nCount ik > nM ax Repetition, exit and apply 
classical algorithm to obtain DCT coefficients; 

If IsSolution = FALSE, exit and apply classical algorithm to obtain DCT 

coefficients; 
} 

Example 14 We choose the first 8 number gray value of the first line of image in fig. 2 as 
vector f . That is, 



f — (foo, foi, fo2, ■ ■ ■ , foiY — 
(156, 159, 158, 155, 158, 156, 159, 158) T 

c = (c , ci, c 2 , • • • , c 7 ) = DCT(f) 

((co) 2 ,(c 1 ) 2 ,(c 2 ) 2 ,---,(c 7 ) 2 ) = 
(1.9814e + 005, 0.51531, 1.5063, 0.95824, 
3.125, 2.5846, 2.7437, 4.4418) T 

Let e = 2.0 x 10" 5 

The computation procedure of 1-D QDCT is listed below: 

Step A(Figl7). Input: AE = \ \f\\ 2 = 198151, (3 = AE, a = ^ = 24768.875; 
Output: The set of solutions S = {(k,Ck) \ ct < (c^) 2 < f3} = {(0,c )}. The output is 
the unique solution (0, Co). 



Insert Fig 17 Here 



Step B (Figl8). Input: AE = AE - (c ) 2 = 11, -4?- = 5.55132 x 10" 5 , (3 = AE, 

ll/ll 2 



a = AE = 1.57143; 

Output: The set of solutions S — {(4, C4), (5, C5), (6, cq), (7, c 7 )}. 

We will obtain one of solution with equal probability. Suppose the output is (6, c 6 ). 

Insert Figl8 Here 

Step C. Input: AE = AE - (c 6 ) 2 = 8.2563, 4^ = 4.16667 x 10" 5 , (3 = AE, a = 

ll/ll 2 

= 1.37605; 

Output: The set of solutions S — {(2, c 2 ), (4, c 4 ), (5, c 5 ), (7, c 7 ), (6, c 6 )}. 
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We will obtain a new solution with probability p = |. That is, repeating subroutine 1 
["-] times approximately, we will obtain a new solution. Suppose the output is (7,07). 

Because AE Z (cr)2 = 1.92505 x 1(T 5 < e, stop. 
Il/ll 2 

E. Time Complexity of Quantum 1-D DCT of Image Compression (QDCT1) 

Suppose the number of retained DCT coefficients is M. 

For typical images, many of the DCT coefficients have values close to zero; these coef- 
ficients can be discarded without seriously affecting the quality of the reconstructed image 
(section I-B.3., Table 1). 

The above example shows that the only one retained DCT coefficientthe is Cq if e = 

5.56 x 10 -5 . That is, Co keeps almost information of vector /. Thus, M = 1 (i.e., M <C N) 
By equation (11), the stepO has complexity 

Time(StepO) = t I + t D < 7t M 

The stepl has complexity 

Time(Stepl) = Time(subroutinel) 
The step2 has complexity 

Time(Step2) = tu + to + St a + Stc < 5£a/ 
Thus, the time complexity of QDCT1 is 

Time(QDCTl) = Time(StepO) + 

M{t D + Time(Stepl) + Time(Step2)) 

< Itu + M{7tM + Time(subroutinel)) 

Because M N, by equation (28), we have 

Time(QDCTl) = 0(VN)t M (29) 

V. QUANTUM 2-D DCT 

By formula (10), the 2-D DCT can be expressed as 

C = DFD = (DF)D = GD (30) 

,where 

G = DF = (gij)NxN = (Di • fj)NxN = 

/ B„ \ 

\ fo fl ••• f N-l J 



(31) 
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A. The Data Structure of Computing Matrix G 



2-D DCT is closely related to the discrete Fourier transform. It is a separable linear 
transformation, that is, the result of 2-D DCT may be obtained by first taking transforms 
along the columns of F and then along the rows of that result 0,0]. 

Clearly, G is the result of 1-D DCT by taking transforms along the columns of F. 

All of vectors A and /^are stored in two different tables of database, where < i, j < N. 

An arbitrary record A or /jCan be fetched into registers of CPU according to indices 
i,j by performing load operation U' L . 
Load operation U' L is defined as 



l«) I/?) 1^10)10)10)10) -4 \ a )\0) \m l/i>|A)|0)|0) (32) 

We design oracles B inner and B' inner to compute the inner product between two vectors 
A and fj 



|a)|/3)|i)|i)|^)|A)|0)|0) B — 



l«>li9>|i>|j>l/i>IA>|0>|(A«/i) 2 > (33) 



l«)l/3)b')b'>i2>iA>io)|o) 



Define a function /' as following 



f(i,j) = \nf<*<pi'f i r<p 

I otherwise 

We design other oracle denoted by 0' f that is used to mark the retained elements of G: 



l«)l/3)K)|j)lS)|A)|0)|(A./;.) 2 )3 



) / ' w |«)l/3)l^)U>l/i)IA)|0)|(A./i) 2 ) (34) 



Let (Bi nner ) 1 denote the inverse transform of Bi nne r- 
Let (f/£) _1 denote inverse transform of U' L . 

Definition 15 By equations (32), (33) and (34), quantum 2-D DCT iteration G' DCT is 
defined as 

G'DCT = (2|e')(e'l - I)(U , L )- 1 (B m ner)- 1 , f B m nerU , L (35) 

N-l N-l 

, where] ? >= £ E E KH^- 
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B. Subroutine 2: Finding an Element gi Q j Q G {gij \ a < (gij) 2 < (3} 
Subroutine 2 is similar to Subroutine 1. 
Algorithm 16 Subroutine 2 
Input parameters: a, f3 

Outputs: subscript (io,jo) and corresponding coefficient g io j ; a global Boolean variable 
IsSolution. 

Stepl. Initialize m = 1 and A = |; IsSolution = FALSE; 

Step2. Choose j uniformly at random among the nonnegative integers not bigger than 

m. 

Step3. Generate the initial state 



1 N-lN-l 

^) = ivEEi«)i' 9 )i!)^i )i )i )i ) 



i=0 j=0 

This can be achieved in 0(2 log 2 N) steps using a 2[log 2 N] — bit Hadamard transforma- 
tion. 

Step4. Apply j iterations of G' DCT starting from the state \^o). i.e., 

= (G' DCT y\v ) 

Step5. Perform computation 

\*end) = B' inner U' L \*i) 

Step6. Observe the 3rd, 4th and 7th register: let (io,jo) and coefficient g io j be the 
outcome respectively. 

Step7. If g io j G {g^ \ a < (gij) 2 < (3}, let IsSolution = TRUE, and exit; otherwise, go 
to step 8. 

Step8. Set m to min{Am, iV} and go to step 2. 

Similar to equation (28), Subroutine 2 has time complexity 

N 

Time(Subroutine2) = 0(— )t M (36) 

yt 

where t is the numberf of elements in set {g^ \ a < (gij) 2 < 13}. 

C. Quantum 2-D DCT (QDCT2) 

Algorithm 17 Quantum Algorithm to compute matrix G approximately 

Input: The line vectors of matrix of D, that are stored in database with indices i (Figure 
10); 

The column vectors of matrix of F, that are stored in database with indices j (Figure 
10); 
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Parameters: AE, a, (3, threshold value e > 0. 

nCounti (0 < i < N): Counting the repeated times of coefficient q that is found when 
repeatedly applying subroutine 2. 

nM ax Repetition: The maximum number of allowed repeated times of subroutine 2. 
nS: The number of solutions that have been obtained. 

Output: the big elements of G and corresponding indices (ia, jo, 9i j ), Ji, <7ijji), • • • , 
(*Kj 3Ri 9i R j R ) sucn that ^ 10 ^ "|j F || 2 ( ' gl - R -'- R ^ < e. The other DCT coefficients are set to zero 
(Table 1). 

N-1N-1 

StepO: AE = \\F\\ 2 = £ E (Afc) 2 ; <* = P = AE '> nS = °! nCounti = 0, where 
< i < A. 

By equation (1) and (2), | \F\ | 2 = 1 1/ 1 1 2 + 1 1 f\ | | 2 + - ■ ■ + 1 |/jv-i| | 2 - That is, we can calculate 
||/o|| 2 ) 1 1 All 2 ; •••! ||/W-i|| 2 one by one , then sum up these values. Thus, 

Time(StepO) = JVtj + (N - l)t A + t D < (5N + 2)t M + (N - l)t A 

while(^ > e) 
{ 

Step 1: Apply subroutine 2 to find a coefficient (i k , j k , g ik j k ); 
Step 2: If IsSolution = TRUE and nCounti = 0, 

then AE = AE - (g ikjk ) 2 , nS = nS + 1; a = /3 = AE; nCount ik = 

nCounti k + 1; 

If IsSolution = TRUE, and nCounti k j k < nM ax Repetition, continue; 
If IsSolution = TRUE, and nCounti k j k > nM ax Repetition, exit and apply 
classical algorithm to obtain DCT coefficients.; 

If IsSolution = FALSE, exit and apply classical algorithm to obtain DCT 

coefficients.; 
} 

Algorithm 18 Quantum 2-D DCT (QDCT2) Appling the above algorithm, we will get 
the result of matrix G. Then compute 2-D DCT coefficients according to equation C = GD 
(i.e., equation (30)) by using the same method. That is, applying the above method two 
times, we will obtain the 2-D DCT coefficients of image. 

We can use classical method to compute the inverse DCT because many DCT coefficients 
is equal to zero so that many running times will be saved. 

Similar to equation (29), Quantum 2-D DCT (QDCT2) has time complexity: 

Time(QDCT2) = 0(N)t M (37) 

In contrast, the classical 2-D DCT has time complexity 0(A 2 log 2 A) multiplica- 
tions. Even the efficiency classical parallel 2-D DCT has time complexity 0(Alog 2 A) + 
Time(Communication) , where Time(Communication) denotes the communication time 
between processors |l7(. Because many data will be shared by all processors, the commu- 
nication time will drastically reduce the efficiency of parallel algorithm in fact By 
contrast, the there is no any communication time will be cost in Quantum 2-D DCT. 
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VI. THE COMPARISON BETWEEN THE QUANTUM DCT ITERATION AND 
THE GROVER ITERATION 



The quantum DCT iteration Gdct = (2|£)(£| — l){Uh)^ x {Pi nner )~ x O ^Oi nner UL is derived 
from Grover iteration G = (2\^>)(^>\ — 1)0. It includes two oracles; by contrast Grover' 
has only one oracle. In classical computer, complex computation must be decomposed into 
many simple function blocks (i.e., oracles) and then puts the function blocks back together 
to construct the complex computation. This is a necessary processing because complexity 
system is always constructed by many simple components The law is also suitable to quantum 
computer maybe. Quantum DCT algorithm in this paper can do two things, which are 
simultaneously calculating DCT coefficients and simultaneously marking the wanted DCT 
coefficients. That is, complex computation has been decomposed into two oracles O inner and 
Of; by contrast Grover's algorithm only can search a record with a given index. Equation 
(25) is the key of applying two or more oracles to perform complex computation. If there is 
no equation (25), applying two or more oracles is no significant (Figl3) 

The method presented in this paper generalizes Grover's algorithm. The method is 
suitable to two or more arbitrary oracles. Clearly, if only if the iteration includes the 
unitary with the form (• • • ) _1 0/(- • • ) such that state is eigenstate with eigenvalue (— 1)^, 
the method is valid. 

The following table shows the difference between the quantum DCT iteration and the 
Grover iteration. 

Insert Table 2 Here 



VII. CONCLUSION 

Discrete Cosine Transform (DCT) is very important in image compression. Almost all of 
digital films such as VCD and digital pictures such as JPEG files are compressed by utilizing 
DCT currently. 

Real-time compressing and high compression ratio are the main research topics of image 
compression. The classical 1-D DCT has complexity 0(N\og 2 N) for iV — dimensional 
vector. The time complexity of the classical 2-D DCT is 0(N 2 \og 2 N) for N x N image. 
In general, the more N is large, the more the compression ratio is high if the whole input 
data sequence is Markov chain and the order of chain is bigger than N. However, when N 
is large, classical DCT does not satisfy the requirement of real-time compressing. That is, 
finding a fast algorithm for large N is significant. 

Quantum Fourier transform (QFT) and Klappenecker's Quantum DCT cannnot be ap- 
plied to image compression directly. 

This paper presents the quantum DCT iteration Gdct that is defined as Gdct = 
(2|£)(£| — I)(L r L) _1 (C hnner)~' l O fO inner U L . And constructs 1-D DCT based on the iteration 
with time complexity 0(y/N) for N — dimensional vector and quantum 2-D DCT with time 
complexity O(N) for N xN image. Two properties of DCT is utilized to design the quantum 
DCT algorithm in this paper. One property is that DCT is energy conservation transform. 
The other property is that many of the DCT coefficients have values close to zero; these 
coefficients can be discarded without seriously affecting the quality of the reconstructed 
image. 
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In classical computer, complex computation must be decomposed into many simple func- 
tion blocks (i.e., two or more oracles) and then puts the function blocks back together to 
construct the complex computation. The law is also suitable to quantum computer maybe. 
The quantum DCT iteration Goer includes two oracles so that it can perform more com- 
plex search; by contrast the Grover iteration G = (2\^f)(^\ — 1)0 has only one oracle so 
that it can only perform simple search. Most quantum algorithms have only one oracle cur- 
rently, such as Shor's algorithm, Grover's algorithm, Simon algorithm, and Deutsch-Jozsa 
algorithm. 

Gdct is a rotation acting on subspace (Fig. 14 and Figl5); by contrast, the Grover itera- 
tion acts on global space. Acting on subspace reduce search range drastically. 

The method presented in this paper generalizes Grover's algorithm. Clearly, if only if the 
iteration includes the unitary with the form (■ • • )~ 1 Of(- ■ ■ ) such that state is eigenstate with 
eigenvalue (— 1)-^ (equation (25), Figl3), the method is valid. The method is universal and 
can be applied to discrete Fourier transform and other transforms for image compression. 
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VIII. LIST OF TABLES 

Table 1: The 2-D DCT coefficients of left-top 8x8 blocks in image (a): Discards all but 
10 of the 64 DCT coefficients in each block, and then reconstructs the image using the 2-D 
inverse DCT of each block. Although there is some loss of quality in the reconstructed image, 
it is clearly recognizable, even though almost 85% of the DCT coefficients were discarded. 
This property is utilized to design quantum algorithm in this paper 



4.9211 -0.00773 0.002581 0.001847 ■ • 





0.014981 -0.00281 0.002071 





0.008015 0.002014 








-0.0196 


















Table 2: The Comparison between quantum DCT iteration and the Grover iteration 
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167 167 168 ■■■ 
165 165 166 ■■■ 

55 60 40 ■■- 

FIG. 1: The image with size 128 x 128 is represented by the right 128 x 128 matrix . Each element 
of this matrix is called pixel. The pixel value fij is a gray-level value. The value ffj is proportional 
to brightness or energy The more fij is large, the more corresponding point in image is bright. 




FIG. 2: 256 x 256 image. 



Gdct = (2|f)(f I - I)(O tnner U L )- L O f (O inner U L )- 
G= (2|fr)(fr| -1)0 

(JJL)-\O mner y x O s O mner XJ L \a) \(3) \i)\f) |0> |0> |0> 

= (-l)/W|a>|^>|i>|7>|0>|0>|0>; 

Q\i) = (-l) m \i) 

G dct has two oracles Oi nner and O j ; 

G has only one oracle 

Gdct is a rotation on subspace; 

G is a rotation on global space. 

Gdct can perform more complex search; 

G only can perform simple search 

The method presented in this paper generalizes the Grover iteration. 



IX. LIST OF FIGURES 
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FIG. 3: reconstructed image 
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FIG. 4: The Parallel Circuit for Computing Inner Product: when N < 2 AX2b , t! < 5t M 
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0i^|O> QFT\\) QFT±=<\0) + \1)) 



FIG. 6: QFT only acts on basis states |0) and |1): QFT-LflO) + |1)) = -h={QFT\Q) + QFT\1)) 




FIG. 7: The schematic diagram of rotating the state \^f) by 6 radians per application of G 



Logical 
Address 



Memory 






1 




N-l 




1 1 

Mapping 


r 


*[0] 


m 




x[N-l] 



FIG. 8: The Conception of the Logical Mapping. The mapping associates data with the corre- 
sponding logical address 
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FIG. 9: The Illustration of the Physical Realization of the Logical Mapping Accessing data 
is very very fast operation so that the time of access can be ignored when designing algorithm. 
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FIG. 10: The Representation of Image by Using Quantum States: LOAD operation Ul '■ \i)\0) 

© fi) is a CNOT operation and is a very very fast operation so that the time of addressing 
can be ignored when designing quantum algorithm such as Grover's algorithm. It is clear that the 
most efficient possible algorithm is in this model of computation 
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FIG. 11: The Conception of Oracle 
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FIG. 12: The Conception of Oracle Of 
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FIG. 13: (Ui) l {Oi n ner) 1 OfOi nner llL only changes the phase of the state as (— This is a 
key of applying two or more oracles to perform complex search. 
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FIG. 14: The Action of A Single Quantum 1-D DCT Iteration Gdct- The initial state |* ) is 
rotated in the subspace span{\i)} by 9 towards the superposition |t) of all solutions to the search. 
Initially, it is inclined at angle | from \a), a state orthogonal to |r). The product operation 
(Ul) -1 (Oinner) -1 fOi nner llL reflects the state about the state \a), then the operation 2|£)(£| — I 
reflect it about |£) in the subspace span{\i)}. 
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FIG. 15: Gdct is equivalent to a rotation on subspace span{\i}} 
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FIG. 16: Schematic processing of subroutine 1 
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FIG. 17: Solution Set: 5 = {c }.The output is (0,c ). 
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FIG. 18: Solution Set S = {(4, C4), (5, C5), (6, (7, C7)}.We will obtain one of solution with equal 
probability.The output is (6, cq) 
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